Non-Fermi liquid behavior of quasi one-dimensional pseudogap materials 
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We study the spectral function of the pseudogap phase of quasi one-dimensional charge density 
wave materials. Using a stochastic approach and emphasizing an exact treatment of non-Gaussian 
order parameter fluctuations we will go beyond a usual perturbative calculation. Our results give a 
good fit to ARPES data and explain the absence of the Fermi edge in charge density wave materials 
even above the Peierls transition, indicating non-Fermi liquid behavior. 
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Non-Fermi liquid behavior is a central topic of mod- 
ern condensed matter physics. A paradigmatic example 
represents the universality class of Luttinger liquids [ij, 3| 
whose spectral properties have been calculated exactly j] 
and which applies to the low-energy physics of a variety 
of one-dimensional (Id) gapless systems. However, re- 
cent angle-resolved photoemission spectroscopy (ARPES) 
measurements of charge density wave materials like the 
blue bronze K0.3M0O3 or (TaSe4)2l H iS,, 2, show 
that these materials cannot be understood in terms of 
Luttinger liquid theory, even at temperatures far above 
the Peierls transition. It turns out that the pseudogap 
observed in the charge density wave materials is similar 
to that found in the underdoped phase of high temper- 
ature superconductors P, ^| . It is this analogy which 
makes a better understanding of the Id pseudogap phase 
even more interesting. 

The relative proximity of the pseudogap phase of 
charge density wave materials to the Peierls transition 
indicates that strong electron-phonon interactions should 
play a dominant role. In the early seventies, Lee, Rice, 
and Anderson studied the single-particle Green func- 
tion of a Id model based on order parameter fluctua- 
tions within the second-order Born approximation. Al- 
though this calculation could qualitatively explain the 
observed pseudogap phenomenon j its validity hinges 
on its perturbative nature. As is well known for corre- 
lated Id electronic systems, interference effects due to 
perfect nesting between the two Fermi points invalidate 
a perturbative calculation near the Fermi level. There- 
fore, one should rather try to find a more accurate or 
even exact solution. A widely used partial fraction ex- 
pansion of the single-particle Green function presented 
by Sadovskii was intended to be exact for Gaussian 
order parameter fluctuations but turned out to be only 
an approximate solution to the problem raised by Lee, 
Rice, and Anderson 

Numerical simulations based on exact diagonalization 
methods 0| or the phase formalism flil, Il7| gave exact 
results for the density of states (DOS) for Gaussian or- 
der parameter fluctuations but showed also that Gaussian 
order parameter fluctuations are inapplicable at low tem- 



peratures where amplitude fluctuations should get frozen 
out and a hard gap should appear in the limit T — > if 
quantum fluctuations are neglected. As pointed out in 
Ref. IitI the right low-temperature physics can be de- 
scribed within a model based on classical phase fluctua- 
tions. Using the phase formalism outlined in Ref. [H the 
DOS was also calculated numerically for the crossover 
regime to non-Gaussian order parameter fluctuations us- 
ing a standard Ginzburg-Landau theory :18]. How to 
flnd non-perturbative results for the complete spectral 
function which could be compared with ARPES data, 
however, remained an open question. 

In this work we shall turn to this problem and calculate 
the spectral function of the Id pseudogap phase of charge 
density wave materials for order parameter fluctuations 
given by a generalized Ginzburg-Landau functional. We 
generalize a stochastic method originated by Halperin 
JMi, 2Ql for a particle in a white noise potential towards 
pseudogap systems involving disorder of the colored noise 
type. 

As a starting point, we use the fluctuating gap model 
[ill 0, 1131 which couples the electronic degrees of 
freedom to a complex order parameter fleld A(a;) = 
A(a;)e'''(^) with amphtudc A(a;) and phase i9(x). The 
two-component wave functions ^pn (x) satisfy a Dirac-like 
equation [2^. 

[-icTsdx + A{x)cr+ + A*{x)(7^]i)n{x) = LL!ni>n{x) ■ (1) 

Here cr± = {exi ± zcr2)/2 and 0-3 are the usual Pauli ma- 
trices. Since the Hamiltonian corresponding to Eq. 
is form invariant under the charge conjugation operation 
H criH*(Ti it is possible to use boundary conditions 
which imply that all eigenstates tpnix) are invariant un- 
der the charge conjugation operation i/'„(x) criip'^{x). 
All these eigenstates may be found by making the Eikonal 
ansatz 



■ip{x) 



_jg-«v(3:)/2 



(2) 



In the thermodynamic limit i — > 00, where L is the 
length of the chain, all physical quantities we will be in- 
terested in are independent of these boundary conditions 
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[20|. The logarithmic amphtude (^{x) and the phase (p{x) 
satisfy the foUowing equations of motion [l7| , 

dxf{x) = 2iLj + 2A(x) cos [ip{x) — i^ix)] , (3a) 
dxC{x) = A{x) sin [ip{x) - . (3b) 

(p{x) and C(a;) are simply related to the DOS p{uj) and 
the inverse localization length (or Lyapunov exponent) 
e-^iio) by [13 

piu;) = ^d^idM^)) , e-\cj) = (a,C(a^)) • (4) 
zvr 

Here and in the following, (...) indicates averaging with 
respect to the probability distribution of A{x) which we 
will specify below. 

The spectral function for right- or left-moving elec- 
trons {a = ±1) is defined in terms of the appropriately 
normalized wave functions ■0„(x) as 



A°'{k;uj)=A{a{kF + k);ij) 

f.L/2 




dx 



L/2 



6{UJ - LJ„) 



(5) 



Note that we measure the deviation fc of a wave vector 
p — a{kp + k) from one of the two Fermi wave vectors 
±kp locally outwards such that due to inversion symme- 
try A°'{k]ijj) is independent of a. 

Instead of using J \tpn\^ = 1 to normalize the wave 
functions let us now use the normalization condition 
|V'n(0)| — 1. Since the process of averaging restores trans- 
lational invariance, it follows 



A%k-u)^{Y. 



-L/2 



5{UJ - Un) 



(6) 

Of course, not every solution to Eqs. ||2Jl and Q is an 
eigenstate of the system. To construct the eigenstates, 
we follow Ref. and integrate the equations of motion 
from the boundaries at —L/2 and L/2 to x — where 
we try to match these solutions. Let us call the solutions 
to Eqs. H3a|) and (|3b|) when advancing the solution from 
=FL/2 to the right (left) '^±{x) and C,±{x)- i^ix) is then 
given by 



tpix) 



je-C+(o)+C+(^)+*v+(^)/2 ^ ^ < 

^^-C-{0)+C-ix)+iv>^{x)/2 ^ 2- > 



(7) 



The eigenstates tpn{x) are determined by the energies 
LOn for which the matching condition <(5+(0) — ^-{Q) (up 
to multiples of 27r) is satisfied. In this case ^p+{x) = 
f-{x) for all X, and the eigenstates have their maxima 
somewhere in the bulk from where the wave functions fall 
off exponentially. 

Let us also introduce 



U±ix) 



ikx px 

dx' e-''="V(x') , (8) 



-0(x) 



which satisfy the additional equation of motion 

d^U{x) ^i[k-u;- A(x)e-*['^(^)-''(^)l]C/(x) + 1 . (9) 

It is now possible to rewrite the spectral fimction as 

^"(fc;c.) = ([[/+(0) - C/_(0)] d{ip+iO) - ^_(0))) 

= 2Re(C/+(0)(5(^+(0)-^_(0))) . (10) 

One may explicitly verify that the spectral function 
A{k;LLj) — A"{\k\ — kp;uj) satisfies the usual sum rules, 



dujA{k;uj) = l, (11) 

rik 

— A{k; cj) = 2(5(^+(0) - ^_(0))) = p{u;) . (12) 



The order parameter fluctuations of A (x) = A(x)e*''^^'' 
are determined by a free energy functional which has 
been derived in the context of superconductivity ''2 1*1 . Ex- 
panding the free energy functional in powers of the order 
parameter field A(x) and derivatives thereof one obtains 
the usual Ginzburg-Landau functional. A better approx- 
imation which is not restricted to small A and hence is 
valid at arbitrary temperatures can be obtained by ex- 
panding only in gradients of the order parameter field. 
The generalized Ginzburg-Landau free energy functional 
including all terms up to second order in the gradients of 
the phase and the amplitude of the order parameter has 
the form El 



rL/2 

/3F{A} = dx 

J -L/2 

, n(A(x)) 



{d^^{x)y + v{Aix)) 



where 
M(A) = 



ii„>0 



ii„>0 



((^2 + A2)5/2 ' 

A2 

((^2 + A2)3/2 ' 



0<cD„<£o 



ttAT 



(13) 

(14a) 
(14b) 
(14c) 



'tL/2 



Here, T = 1//3 is the temperature of the system, A <C 1 
is the dimensionless electron-phonon coupling constant, 
eo is an ultraviolet (bandwidth) cutoff, and On = (2n -I- 
1)t:T are fermionic Matsubara frequencies. In the weak 
couphng hmit A <g; 1, V{A) and hence i^{A} depend 
on A and eo only via the critical mean-field temperature 
T^F = 1.134 eoexp(-l/A). 

The probability for a given realization of the order 
parameter field A(x) is proportional to exp(— /3F{A}). 
This probability can be described as a stochastic pro- 
cess which due to the locality of F{A} has the Markov 
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FIG. 1: Temperature dependence of ^J{K^xy) and ^ ^ for a 
lattice spacing d = 0.2 vf/T^^ such that C ^ d for T < T^^ . 



property. The probability P(A) that A{x) assumes the 
value A and the transition probabilities T{A{x) ^A jx') ) 
can be calculated using the transfer matrix method |22l | . 
At this point we would like to note that the integral 
and derivatives in Eq. H13|l should be considered as a 
sum and finite differences. In the case of a generalized 
Ginzburg-Landau theory as considered here, P{A) and 
T{A{x) ^ A{x')) actually do depend on the lattice con- 
stant d. In Fig. n we show ^ (A^) and as functions of 
temperature. For small temperatures we have ^ ~ vp/ttT 
and (A^) w 1.76 T^^, consistent with the predictions 
of a strictly Id theory based on classical phase fluctua- 
tions miijl. Since ^(A^) > 1.76 T^*^^, our generalized 
Ginzburg-Landau functional given in Eqs. I|13|) and l|14l) 
may not simply be replaced by a usual Ginzburg-Landau 
functional by expanding in A. In contrast to predictions 
based on standard Ginzburg-Landau theory as used in 
Refs.ll^andE3. our calculations show that (A^) grows 
with increasing temperature. 

The stochastic process corresponding to T{A{x) <— 
A{x')) can be described by the Langevin equations 



d^A{x) ^ a{A) + b{A) r]A{x) 



d.^ix) = 



n{A{x)) 



1/2 



ri^{x) 



(15a) 
(15b) 



Here 77^(0;) with i = A^-d represent Gaussian white noise 
with {r],{x)r]j{x')) S,jS{x - x'). a(A), 5(A) and P(A) 
are simply related to r(A ^ A') and satisfy a(A) = 
6(A)9a(KA)P(A))/2P(A). 

To make further progress we absorb the phase of the 
order parameter field i^{x) by the gauge transformation 
IIIII2JI ^(x) (fix) + i?(a;). We now obtain the equa- 
tions of motion 

dxf{x) — 2uj — dx^{x) + 2A{x) cos [^{x)] , (16a) 
dxC{x)^ A{x)sm[^{x)] , (16b) 
dxU{x) ^i[k-LU- A(x)e~''^(=")]f/(x) + 1 . (16c) 

Since the stochastic process for a given realization of the 
order parameter field has the Markov property we may 



write Eq. lfTU|l as 

A"{k;uj) =2Re / dA [ 
Jo Jo 



dip ■ 



1 



P(A) 

xPi(A,(^;0)Po(A,-(^;0) . 



(17) 



Here 

Po{A,ip;x) 
Pi(A,(^;x) 



{Siifi - ip+ix))S{A ~ Aix))) , (18a) 
{U+{x) 5{ip - (p+{x)S{A - A(a;))} . (18b) 



An expression similar to Eq. H17|l (with Pi replaced by 
Pq) does hold for the DOS. Making use of the Gaussian 
white noise property of the variables rji (x) , one may show 
that Pq and Pi obey (with Oa acting on everything to its 
right) 



dxPo = ^ 2(9^(cj + A cos c/j) , 2 



- ^AajPo , 
(19a) 



9.P1 



92 

{dAbf 



2{uj + A cos ip)d^ +i{k-u}- Ae**^) 
aAa]Pi+Po. (19b) 



The stationary solutions to the 2d Fokker-Planck equa- 
tion H19a|) and the 2d partial differential equation (jl9b|) 
may be found numerically. We then obtain the spectral 
function A°'(k;u!) (as well as the DOS and the inverse 
localization length) by simple quadrature. The spectral 
function, the DOS and the inverse localization length are 
therefore simultaneously directly calculable without any 
approximations. There is no need to simulate a long 
chain as was done for the DOS (and the inverse localiza- 
tion length) in Refs.ElEEllii 

As can be seen in Fig.|21 at low temperatures A°'{k; lu) 
exhibits two peaks instead of a quasi particle peak. 
Although the qualitative features of the spectral func- 
tion agree wi th p erturbation theory or the Sadovskii ap- 
proximation |25| . on a quantitative level there are large 
deviations, being due to the non-Gaussian nature of the 




FIG. 2: Spectral function A°'{k;uj) for fc = and k = 
0.5T^^/vf. The temperature varies form T/T^^ = 0.1 to 
1.0 in steps of 0.1. 
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FIG. 3: (a) A<{k;io) for T = 0.5T^^ and vpk/T^^ = -2 
to 2 in steps of 0.25. Thick line: A: = 0, i.e. p = kp. The 
spectra are normalized to the same maximum hight. (b) Solid 
lines: p^{uj) = / ^A<{k;co) for T/T^^ = 0.3 to 0.8 in steps 
of 0.1. Comparison with experimental data: Empty circles 
(Ref. 0): T = 250K._Empty triangles (Ref. 6): T = 300 K. 
Filled triangles (Ref. 4): T — 313 K. As a fit parameter we 
have used T^^^ = 60 meV= 696 K. 



order parameter fluctuations. To compare our results 
with experiments, we note that for large photon energies 
the conventional interpretation of the ARPES line shape 
predicts that it is essentially proportional to 0, 



(20) 



where A^{k; uj) = A°'{k; ui) / (e^'^ + 1) is the spectral func- 
tion multiplied with the Fermi function. In Fig. ISJa) we 
show a plot of A<{k;uj) for T = 0.5T^^ and different 
values of k. A^{k;uj) agrees reasonably well with the 
ARPES spectra presented in Fig. 1(a) of Ref. which 
were taken at T = 210 K. In Fig.|3Jb) we compare angle 
integrated photoemission data with our theory. The only 
fit parameter used is = 60meV= 696 K. This value 
is less than a factor of two larger than « 400 K, as 
obtained from optical conductivity measurements |23l |. 
The deviation might be due to a renormalization of A by 
electron-electron interactions and the simplicity of the 
model. Nonetheless, {lj) — J {k; lj) does give a 

good fit to experimental data. We conclude that order 
parameter fluctuations seem to dominate the pseudogap 
phase. 

In summary, we have calculated the spectral function 
of a Id electronic system coupled to order parameter fluc- 
tuations determined by a generalized Ginzburg-Landau 
theory. Our results show that the essential physics of the 
pseudogap phase of quasi Id charge density wave systems 
can be captured within a strictly Id model based on or- 
der parameter fluctuations, treating their non-Gaussian 
statistics exactly. These order parameter fluctuations 



lead to non- Fermi liquid behavior and the experimentally 
observed absence of a Fermi edge. 
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